0^ 



Capture on High Curvature Region: Aggregation of Colloidal 
Particle Bound to Giant Phospholipid Vesicles 

Quan-Hui Liu^'^H, Ji-Xing Liu\ Ou-Yang Zhong-Can^ 
^Institute of Theoretical Physics, Academia Sinica, P.O. Box 2735, Beijing, 100080, China 



^ \ Department of Physics, Hunan University, Changsha, 410082, China 

O ; (February 1, 2008) 

o; 

^ . Abstract 

-)— » . 
c^ . 

S: 

I ' A very recent observation on the membrane mediated attraction and or- 

^ I dered aggregation of colloidal particles bound to giant phospholipid vesicles 

I ^; (I. Koltover, J. O. Radler, C. R. Safinya, Phys. Rev. Lett. 82, 1991(1999)) 

►^ I is investigated theoretically within the frame of Helfrich curvature elasticity 

(n' 

l/^ ' theory of lipid bilayer fluid membrane. Since the concave or waist regions of 

^ I the vesicle possess the highest local bending energy density, the aggregation 

G^ ■ 

Q>^ I of colloidal beads on these places can reduce the elastic energy in maximum. 
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C^ ' Our calculation shows that a bead in the concave region lowers its energy 

B' 

I ' ~ 20kBT. For an axisymmetrical dumbbell vesicle, the local curvature en- 

TD ■ 
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Q I ergy density along the waist is equally of maximum, the beads can thus be 
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^ I distributed freely with varying separation distance. 
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Recently, much attention focused on the experimental results carried out by a group 
in University of California (GUOC) on the interaction between DNA, colloidal particles 
and lipid membrane |l|-§[. These experimental advances are helpful to understand the 
structure of bioassembles, the possible vesicle shape transformations, the capture function 
and the transport processes of the membrane, etc.; and all of these have suggested their 
potential use for industrial as well as medical applications |Q. In 1998, even two independent 
theoretical groups gave the same theory to explain the structure of DNA-cationic lipid 
complexes observed by the GUOC [|],^. In March of 1999, the GUOC reported a light 
microscopy study of the interactions between colloidal beads either chemically bound or 
physisorbed onto flexible giant vesicle fluid membranes 0. But the chemically bound cases 
were a continuous work of their previous one |jl|J^ . Their new experiments were mainly on the 
physisorption between the colloidal beads and membranes. They observed two essentially 
new phenomena via two types of different experiments. 

As GUOC reported, the initial bare unilamellar vesicles of sphere topology were prepared 
from phospholipids POPC (l-Palmytoyl-2-Oleoyl Phosphocholine) or DMPC (Dimiristoyl 
Phosphocholine) and 0.5% biotin-X DHPE (molecular probes) following a standard protocol 
101, and the typical vesicle has a radius R ~ IQ^m. The colloidal beads were chemically 

grafted with treptavidin a water soluble protein which specifically binds to biotin; and 

the typical bead has a radius a ~ 0.45/im. After putting the beads to the surfaces of vesicles, 
treptavidin beads were linked to giant vesicle doped with biotinilated lipids by mixing dilute 
suspensions of each component. First type of their experiment is on the shapes of equilibrium 
bead-membrane hybrid vesicles. They found that these vesicles clearly revealed concave 
regions where the beads locate. For further demonstration of this fact, they had captured 
vesicles using suction glass pipettes then released them, and found that the vesicles reverted 
to the nearly original shapes with the bead centered at concave regions again. We noticed 
that when the bead-membrane hybrid vesicles are either in free states or in head of the 

suction glass pipettes, the bead pinched in membrane and the degree of pinch which 

is balanced by the adhesion force does not appear significant change. When there are 



many beads on vesicles, these beads either segregated condensed phases of beads located 
in concave regions of membranes, or firmly pinned to the waist of a nearly axisymmetrical 
dumbbell shape vesicle along the meridional direction but distributed freely with varying 
separation distance along the azimuthal direction. The second type of their experiment is on 
the movement of the two beads having the same radius on a quasi-spherical vesicle. At first, 
these two beads walk randomly in the surface of the vesicle. If their center-center distance 
was within four times of their radius, an attraction drives them together. These two beads 
finally formed a concave region and remained there from then on. 

From their experiments, two conclusions seem to be definite. First is that there may be 
no attraction even on the membrane otherwise the beads in the waist of an axisymmetrical 
dumbbell shape vesicle could be condensed. Of course, there is no evidence to show that there 
are observable repulsion between them. Second, whenever the vesicles adhered by beads are 
in free states or are captured by suction glass pipettes, and whenever the beads are randomly 
walking on of the membranes or trapped in concave region, there is no observable completely 
enclosed or laid bare phenomena for the pinched beads. It means that the adhesion energy 
of each bead is nearly the same as long as it is on the membrane. We have noted that 
the GUOC has proposed a theoretical explanation for such phenomena based on a known 
theory developed by Dan 0, and indeed offered many useful theoretical information. But 
their explanation is not sufficient enough to put all of their results into a theory coherently. 
For example, from the Dan's theory, the membrane-induced attractions are caused by the 
adsorbing DNA perturbing the equilibrium packing of the lipids, and the range of interaction 
is typically of order nm. But here the distance of two beads starts attracting is LSfim and 
each bead at least four orders heavier than the DNA molecula, and there is no attraction 
between beads in the waist line of an axisymmetrical dumbbell. Whether we can apply the 
principle of Dan's theory to our situation seems to be a problem. Moreover, the concave 
regions seem not to be represented by catenoid shapes as did in 0, because the mean 
curvatures and therefore the curvature energies of such shapes are zero, and any further 
deformations of the membranes will increase the elastic energy. Thus the localization of 



beads on concave regions is energetically unfavored. In this Letter, we will present an 
improved theory to account for all these phenomena. In our approach, the concave or waist 
regions of the vesicle possess the highest local bending energy density, the aggregation of 
colloidal beads on these places can reduce the elastic energy in maximum. An explicit 
calculation will show that a bead in the concave region lowers its energy ~ 20kBT where ks 
is the Boltzmann constant and T the temperature. For an axisymmetrical dumbbell vesicle, 
the local curvature energy density along the waist is equally of maximum, the beads can 
thus be distributed freely with varying separation distance. 

First of all, we take for grant, as experiments implied (a quantitative calculation will 
be given latter), that the area of pinching membranes around each bead of the same radius 
is fixed and the adhesion energy for each bead then is regarded as the same as long as 
the beads are on the membrane. The stable shapes of bead-membrane hybrids are of course 
energy minimum states, and the two energy degenerate shapes of sphere topology vesicle are 
axisymmetrical biconcave shape and the axisymmetrical dumbbell shape from the instability 
analysis of a sphere shape vesicle [^]. So, the formations of concave regions or waist for sphere 
can be taken place either before the beads were introduced, or while the bead-membrane 
hybrids look for the energy minimum states. According to experimental facts 0, if there 
only one bead, the bead will locate on the concave region; if there are many beads on it, 
there must be a bead at the center of the concave region first for that a bead-membrane 
hybrid vesicle after all develops a concave region to trap these beads. The main text of this 
Letter is devoted to show why a bead favors the concave regions and how two beads attract 
in detail. For our purpose, we will use an exact solution representing a concave shape in the 
fluid membrane to show that the free energy density per area in the vesicle is not uniformly 
distributed; and seriously deformed areas, concave and waist regions for example, have the 
much higher energy density than elsewhere. Location of the beads in these places can make 
the whole energy lower down. The membrane mediated attraction is thus induced by the 
free energy difference. 

In physics, equilibrium shapes of phospholipid vesicles are assumed to correspond to the 



minimum of the elastic energy of the closed bilayer membrane consisting of the amphiphilic 
molecula. The energy functional or the total free energy F of a vesicle is given by Helfrich 
spontaneous curvature theory as 0: 

F = -k^ f{ci + C2- cofdA + 6p fdV + X fdA, (1) 

where dA and dV are the surface area and the volume elements for the vesicle, respectively, 
kc ~ SOksT is an elastic modulus [P,|^, ci and C2 are the two principal curvatures of the 
surface, and cq is the spontaneous curvature to describe the possible asymmetry in both sides 
of the bilayer membrane. The Lagrangian multipliers 6p and A take account of the constraints 
of constant volume and area, which can be physically understood as the osmotic pressure 
between the ambient and the internal environments, and the surface tension, respectively. 
In GUOC experiment, the tension energy X J dA and the pressure energy 6p J dV appear 
not to change significantly, then we will simply discard these two constant terms hereinafter. 
Here we have not included the interaction between the beads and membranes, but will treat 
them as perturbations. 

Variational calculus of the elastic energy Eq. (|I]) gives the equilibrium shape equation [Q : 

(2H + Co) {2H^ -2K- coH) + 2\/^H = 0, (2) 

where v^ = ~i=^i{9^^ \/9dj) is the Laplace-Beltrami operator, g is the determinant of the 
metric Qij and g'^^ = {gij)~^, K = ciC2 is the Gaussian curvature and H = — (l/2)(ci + C2) is 
the mean curvature. Assuming that the shape has axisymmetry, the general shape equation 
Eq.(2) becomes a third order nonlinear differential equation ||^: 

COS ip{—-ir) = Asm ip COS ip[— —-)[——)— COS ip [sm ip cos ip){—r) 

dr'^ dv^ dr 2 dr 

7 sm-ip cos'^ ifj dij) 2 2 cos^ ^ d'^^jJ 

2r dr r dr"^ 

,c^ 2coSinV' siia^ ih — 2 cos"^ iIj , ..dih, 

+ If -r S^ 1 "= ^(^) 

.c^smih sin^'?/' + 2sin^cos^'?/', . . 

+ l^, ^3 )■ (3) 



where r is the distance from the symmetric z axis of rotation, iplr) is the angle made by 
the surface tangent and the r axis. The positive direction of the angle is that the angle 
measured clockwise from r axis. It is contrary to the usual mathematical convention; so, 
the mean curvature H is — l/2(sin?/'/r + dsimp/dr), in which ci = sinip/r denotes the 
principal curvature along the parallels of latitude, and C2 = dsinip/dr denotes that along 
those of meridians. To note that different papers may use different sign conventions, and 
the sign convention used in this Letter is self-consistent. It is worthy of mentioning that the 
spontaneous curvature Cq carries a sign. When the normal of a surface change its direction, 
ci,C2 and cq must change their signs simultaneously. Keeping the directions of r and z as 
the usual, we have therefore: 



"■^ 



dz/dr = —taiiih(r) 

(4) 

z{r) — -2(0) = — Jq ta.iailj(r)dr. 

We see that the sphere with parameterization simp{r) = v/Rq, where i?o > is the radius 
of the sphere, is always a solution of Eqs. (0,|^). It is easily to verify that Eqs. ([I|,0) have 
the following solutions [|T^] 

sin ip = r + CqV In r, (5) 

where we take r is measured with lOfim as the length unit and so for the value of spontaneous 
curvature Cq. When cq = 0, this solution gives nothing but sphere of radius lO/im. When 
Co is a small quantity, positive cq leads to oblate ellipsoid and negative Cq leads to prolate 
ellipsoid respectively |]T2[. Increasing cq from zero, we can have a great variety of the concave 



regions from shallow to deep. The elastic energy density, the energy per area, from Eqs. 



is: 



dF/dA = 2kc{l + colnrf. (6) 

To obtain a typical concave region, an average shape of concave regions trapping beads 
in the Figs.l, 2, 5, of GUOC report |^ is used for fitting, and such a concave region is 
shown in Fig.l by the upper half of a sphere of radius 10/xm. Its parameterization is 



sinip = r + 2.1rlnr, and we plot the energy density by thin hne in Fig.2. The elastic 
energy in the whole concave region, i.e. in the upper half piece of the surface of quasi-sphere 
shape, can be easily carried out to be 280.61kBT. The contact area a between each bead and 
membrane are reasonably a = irr^ = 7r(0.3)^(/im)^ as shown in with tq = O.Sfim being 
the projected radius 0, and the adhesion energy is comparable with elastic energy in this 
contact area, which is 2kcCr/a'^ ~ lOOksT [§]. This is convincing because it guarantees that 
the beads can never escape from the vesicle due to the thermal agitation. Since the adhesion 
energy is independent from the place where the beads locate, the more serious deformed the 
original membrane is, the more bending energy is released from the whole free energy. The 
original bending energy in the contact area of the membrane then plays the role of trapping 
energy. Since the beads constrained in the membrane, it has only two degrees of freedom. 
The kinetic energy of each bead is therefore 2{l/2kBT) = ksT. If the trapping energy is 
less than this energy, the beads reveal Brownian walk motion. The trapping energy in the 
center of the concave region is AF(0): 

AF(0)= 2k^ / (1 + 2.llnrfdA 

= 2kc / (1 + 2.llnrf/^l - (r + 2.1rlnr)27rrdr 

= 19.8kBT ~ 20kBT. (7) 

The trapping energies in range r > 2fim are given by AF{r) = 2kc{l + 2.1/nr)^cr = 
27r(0.3)2A;,(l + 2.1lnrf < 0.96fcsT, as shown by the bold line of the Fig. 2. There- 
fore in these interval, the Brownian motion would dominate the bead's trajectory. When 
r = 4 * OAdfim = l.Sfim, AF{1.8fim) = LIS^bT, the trapping energy is comparable with 
the thermal fluctuation. Once the bead is within r < l.Sfim, the trapping energy over- 
whelms the thermal fluctuation and would drive the bead move to the center. Once a bead 
moves to the center, its Brownian motion is energetically depressed and the bead can be 
regarded as trapped. This gives a reasonable explanation to the primary experimantal result 
as the tendency of a bead to locate on a concave region as observed in [^ . 

For further demonstration the degree of agreement of our theory with the experimental 



date, we estimate the magnitude of relative velocity v of the motion of two beads from 
relative position r = 1.8/im to r = 0.9/im with one bead already locating at the center of 
the concave region as shown in Fig. 1. We have approximately l/2mf^ ~ Ifc^T— IO/cbT, i.e., 



we have f ~ f o — 3f o with vq = J2kBT /m. By comparing these date with the experimental 
ones, we see that agreement is twofold. First it means that in whole range r < l.Sfim, there 
is an increase of velocity from vq to 3fo; it is true from the observation as shown in Fig. 3(b) 
1^. Second, the average of interaction strength in the range r = 1.8 fim to r = 0.9^m is 
bhsT, which can be regarded as the quantitative statement of "several fc^T" |^. 

For checking the above theoretical consistence of our approach, we return to verify the 
basic assumption about the constancy of the contact area a = 7r(0.3)^(/im)^, i.e., the con- 
stancy of the adhesion energy. Since the elastic energy has the scaler invariance, all sphere 



shapes have the same energy Snkc ~ TSO/c^T, independent from the radii [|TT|. Therefore, 
if a small change of the contact area as a = 7r(0.3 ± 0.1)^(/im)^ is allowed, it would accom- 
pany the large energy change AF = 37.21kBT — 148.95A;bT. So large energy change on the 
membrane must have the observable transport effect of membrane materials, but no such 
processes were reported in the experiment 0. We would also like to compare the elastic 
energy of a perfect sphere shape vesicle with it having a concave region in its half piece as 
shown in Fig. 1, the elastic energy reduces an amount 376.46A;bT — 280.61/^^^ = 95.85kBT. 
Therefore the formation of the concave regions is also energetically favored. 

Till now we have developed a consistent theoretical approach to explain the primary 
experimental results of GUOC [^: The tendency of a bead to be captured on the center 
of a concave region (the region of high curvature region), and the membrane mediated 
attraction due to the elastic energy difference of the membrane. Such principle can be used 
to explain other observations given in GUOC report 0. If there are four, five and six 
such beads, these beads will closely occupy the corners of a square, pentagon and hexagon 
respective, leaving the center unoccupied, since the these beads tends to occupy as large area 
of membrane as possible. Since along the waist of a axisymmetrical dumbbell shape vesicle, 
the energy density is in maximum and equally distributed in the surface; the beads can thus 
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be distributed freely with varying separation distance. These facts are again evident from 
the experimental observations 0. 

In summary, our approach gives a more consistent mechanism of physisorption than the 
others [0]: the maximum reduce of free energy can induce attraction of colloidal beads on the 
appreciated position of the membrane. Undoubtedly, if the beads are not on the membrane, 
there is no such attraction, which is confirmed with the experimental observation too [^. 
This may have significance in cell biology. In fact, in biophysics, we have known that the the 
absorption of protein and ion channels mainly in the concave regions in the cell biomembrane 
such as red blood cell |T^, but these phenomena have not be sufficiently studied before. 
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FIGURES 




FIG. 1. A concave region ( only the upper half part is relevant ) fitted those given by experi- 
mental ones in Figs. 1, 3, 5 of GUOC's paper[3]. The parameterization is sinip{r) = r + 2.1r Inr 
with 10/im as the length unit. 
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FIG. 2. The dependence of energy per area and the trapping energy with r in a concave 

region sin^ = r + 2.1r In r,r G (0, lO/um). Thin hne shows the elastic energy per area using unit 

5000kBT{fim)^'^. Bold line shows the elastic energy in contact area a = 7r(0.3)^(/im)^ (trapping 

energy) using length unit of 10{pm) in r axis. 



12 



